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Abstract 

Conditions  for  the  turn-on  of  short  wavelength  flute  modes  on  plasmas  expanding 
across  an  external  magnetic  field  are  discussed.  The  experimental  data  concerning  scaling 
of  the  instability  onset  as  a  function  of  magnetic  field  are  reviewed,  followed  by  a  similar 
review  of  numerical  simulation  results.  Theoretical  conditions  for  the  appearance  of  the 
fast  growing,  short  wavelength  lower  hybrid  drift  instability  are  derived.  Overall,  the 
experiments  and  simulations  in  which  a  scaling  can  be  deduced  suggest  the  turn-on  of  the 
instability  is  a  strong  function  of  the  magnetic  field.  These  results  are  discussed  in  relation 
to  the  NRL  laser  experiment,  which  in  contrast  exhibits  little  dependence  of  the  instability 
onset  on  the  applied  field. 
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1.  Introduction 


Numerical  experiments  in  which  a  plasma  is  created  and  then  expands  across  an 
externally  imposed  magnetic  field  show  that  the  surface  becomes  rippled,  forming  field 
aligned  (flute)  modes.  Examples  include  gun  generated  plasmas  [Dickinson  et  al.,  1962], 
laser  produced  plasmas  [Okada  et  al.,  1981;  Ripin  et  al.,  1987],  and  gas  releases  in  space 
[Bernhardt  et  al.,  1987].  While  these  plasma  have  widely  differing  parameters  (as  described 
later),  the  fact  that  they  exhibit  similar  behavior  suggests  a  common  instability  mechanism. 
Numerous  theoretical  and  simulation  efforts  have  addressed  the  origin  of  these  surface 
structures.  Akimoto  et  al.  [1988]  summarize  the  recent  experimental  and  theoretical 
work. 

Presently,  there  is  a  general  consensus  concerning  how  the  instability  is  produced.  On 
the  time  and  distance  scales  of  interest,  the  ions  are  unmagnetized  and  expand  outward 
from  the  source  region  unimpeded  by  the  magnetic  field,  but  constrained  by  the  electrons 
that  are  tied  to  the  field,  via  a  radial  electric  field  that  is  set  up.  The  azimuthal  Ex  B  drift 
of  the  electrons  relative  to  the  ions  is  the  free  energy  source  for  generating  the  unstable 
waves.  Theoretically  it  is  easily  shown  [e.g.,  Winske,  1988]  that  the  range  of  instability 
extends  from  high  frequency,  fast  growing  (~  lower  hybrid  frequency),  short  wavelength 
(electron  inertial  length  or  shorter)  modes  to  much  longer  wavelength  (~  ion  inertial  scale), 
slower  growing  (~  ion  cyclotron)  waves.  This  range  of  scales  has  led  to  different  approaches 
io  derive  the  properties  of  the  unstable  waves  and  to  different  names  given  to  the  instability. 
Hassam  and  Huba  [1987,  1988]  start  with  a  fluid  approach  including  the  Hall  term  and 
obtain  the  long  wavelength  limit  of  the  instability,  which  they  call  the  “unmagnetized  ion 
Rayleigh-Taylor  instability”,  while  Winske  [1988,  1989]  and  Galvez  et  al.  [1988]  employ 
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a  kinetic  treatment  and  concentrate  on  the  short  wavelength  version  of  the  instability, 
here  called  the  “driven  lower  hybrid  drift  instability”.  Both  approaches  have  a  serious 
shortcoming,  in  that  neither  gives  the  correct  scaling  of  the  wavelength  of  the  instability. 
The  fluid  approach  includes  the  wavelength  of  the  instability  as  a  free  parameter;  the  kinetic 
method  predicts  a  wavelength  corresponding  to  maximum  growth  that  is  consistent  with 


simulations  and  some  experiments,  but  is  more  than  an  order  of  magnitude  shorter  than 
that  observed  in  other  experiments  [see  Winske,  1989]. 

Leaving  aside  this  unresolved  issue  of  the  wavelength,  we  turn  to  another  question  of 
interest  at  present,  namely  the  scaling  of  the  turn-on  of  the  instability  with  parameters, 
particularly  the  applied  magnetic  field  and  the  ion  mass.  Huba  [1987]  has  shown  that  the 
long  wavelength  unmagnetized  ion  Rayleigh- Taylor  instability  has  an  onset  condition  that 
is  independent  of  the  magnetic  field,  which  is  in  agreement  with  experiments  carried  out 
with  the  NRL  laser.  The  issue  at  hand  is  what  are  the  corresponding  turn-on  criteria  for 
the  driven  lower  hybrid  drift  instability,  how  do  they  scale  with  parameters,  and  what 


do  they  predict  for  the  experiments.  To  facilitate  the  discussion,  the  rest  of  the  report 
is  broken  up  into  four  sections.  Sections  2  and  3  review  the  relevant  experiments  and 
simulations,  respectively,  to  infer  possible  onset  criteria.  Then,  based  in  part  on  these 
results  as  well  as  the  linear  properties  of  the  instability,  turn-on  conditions  for  the  driven 
lower  hybrid  drift  instability  are  derived  in  Section  4.  Generally,  the  four  onset  conditions  — 

Ft 

that  are  considered  all  show  a  rather  strong  dependence  on  the  applied  magnetic  field,  ! 

,  * 

in  contradiction  to  the  laser  experiments  of  Ripin  et  al.  [1987].  These  differences  are 
discussed  further  in  Section  5,  where  possible  tests  using  experiments  and  simulations  fo 


resolve  the  discrepancy  are  proposed. 
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2.  Review  of  Experiments 

The  results  of  four  different  experiments  in  which  surface  waves  are  observed  are 
discussed:  (A)  the  AMPTE  barium  releases  in  the  magnetotail  [Bernhardt  et  al.,  1987], 
(B)  the  plasmoid  experiments  of  Dickinson  et  al.  [1962],  and  laser  generated  plasma 
experiments  of  (C)  Okada  et  al.  [1981]  and  (D)  Ripin  et  al.  [1987]. 

A.  Preliminaries 

Before  proceeding,  we  introduce  a  few  basic  concepts  that  will  be  used  repeatedly 
throughout.  The  experimental  arrangement  we  wish  to  consider  is  shown  schematically  in 
Figure  1.  The  initial  conditions  consist  of  a  plasma  source  containing  N  ions  and  electrons 
with  initial  radius  expanding  radially  outward  with  velocity  VD  into  a  low  density 
background  plasma  (or  vacuum)  across  a  magnetic  field  of  magnitude  B.  The  expansion 
velocity  Vd  is  assumed  to  be  subAlfvenic  and  the  initial  plasma  beta  0  »  1,  so  that 
as  the  plasma  expands,  it  excludes  the  magnetic  field  and  forms  a  diamagnetic  cavity. 
The  expansion  continues  out  to  the  magnetic  confinement  radius,  Rb,  where  the  excluded 
magnetic  field  energy  balances  the  initial  kinetic  energy  of  the  expansion: 

\Nm.vl  =  fly  (1) 

where  V  is  the  final  volume  (V  =  47rJ2g/3  in  3-D,  ir RgL  in  2-D). 

The  above  simple  argument  assumes  that  the  expansion  is  everywhere  perpendicular  to 

B.  While  this  is  a  reasonably  good  approximation  in  3-D,  one  can  include  the  free  expansion 
along  B  as  well  [Gisler  and  Lemons,  1988].  The  convenient  parameter  to  measure  distance 
is  Rb ■  However,  because  Rb  changes  with  magnetic  field  strength,  which  is  one  of  the 
parameters  we  wish  to  vary,  we  will  imagine  an  experiment  with  reference  parameters 
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Figure  l.  Schematic  of  plasma  expansion  experiments:  N  ions  and  electrons  with  initial 
radius  R ,  expanding  into  magnetic  field  B  with  velocity  Vj>;  expansion  stops  with 
formation  of  diamagnetic  cavity  of  radius  Rb<  c 


B  =  B0  with  corresponding  RB  =  R0  and  use  those  parameters  as  a  basis  for  scaling 
studies.  Similarly,  scalings  will  be  done  with  a  reference  ion  mass,  m0,  and  one  could  carry 
out  a  comparable  analysis  for  other  initial  parameters  as  well. 

From  the  same  dynamical  considerations,  one  can  show  that  the  velocity  as  a  function 
of  times  is 

V{R)  =  VD[l  -  {R/RB)a}1/2  (2) 

where  a  =  2  (2-D)  or  3(3-D).  By  integrating  dR/V(R)  from  R=0  to  R  =  RB,  one  obtains 
the  time  corresponding  to  maximum  expansion  of  the  cloud 

tB  =  aRs/Vo  (3) 

where  a=1.4  in  3-D,  7t/2  in  2-D.  A  third  important  parameter  is  pt,  the  directed  ion 
gyroradius 

Pi  =  vD/ni  (4) 

where  Hi  =  eB/m,c  =  ion  gyrofrequency.  Our  previous  study  [Winske,  1989]  has  suggested 
that  Pi/RB  is  the  single  most  important  parameter  to  characterize  the  strength  of  the 
driven  lower  hybrid  drift  instability. 

B.  AMPTE 

The  AMPTE  (Active  Magnetospheric  Particle  Tracers  Explorers)  mission  [Krimigis 
et  al.,  1982]  involved  a  number  of  chemical  releases  in  space,  one  of  which  occurred  in 
March  1985  with  the  release  of  barium  ions  in  the  magnetotail.  In  this  experiment  about 
1025  atoms  were  released  with  initial  velocity  VB  ~  1.3km/sec  into  the  hot  magnetotail 
plasma  and  ambient  magnetic  field  of  about  8nT.  Because  of  the  finite  ionization  time  for 
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the  barium  (about  30  sec)  and  the  expansion  of  the  cloud,  the  actual  barium  density  is  a 
complicated  funftion  of  time.  Nevertheless,  ground  based  observations  showed  a  maximum 
expansion  radius,  R  ~  210km,  consistent  with  that  calculated  from  (1),  Rb  ~  240km, 
taking  into  account  the  free  expansion  along  the  magnetic  field.  Furthermore,  for  these 
parameters  /?,  ~  210km  as  well.  Of  interest  here  are  the  optical  observations  of  Bernhardt 
tt  al.  [1987],  who  saw  striations  on  the  surface  of  the  cloud  starting  at  tc  ^  180sec, 
corresponding  to  Rc  <  Rb,  which  persisted  beyond  the  time  of  maximum  expansion, 
ts  —  240sec,  and  the  beginning  of  the  recollapse  of  the  cloud.  From  the  optical  images, 
it  was  determined  that  the  mode  number  of  the  instability  was  m=24,  which  corresponds 
to  a  wavelength  A  ~  55km,  using  the  measured  R=210km  at  maximum  extension.  The 
wavelength  is  considerable  longer  than  that  corresponding  to  the  driven  lower  hybrid  drift 
instability,  as  discussed  recently  by  Winske  [1989]. 

C.  Plasmoid  Experiments 

The  second  set  of  experiments  to  be  considered  were  carried  out  by  Dickinson  et  al. 
[1962],  in  which  a  plasma  was  produced  by  directing  a  plasma  gun  against  a  copper  barrier. 
The  resulting  plasma  expands  with  velocity  V o  ~  104m/sec  across  a  magnetic  field  that 
could  be  varied  in  strength  from  0  to  0.8T.  In  this  case  the  composition  of  the  plasma  is  not 
known  very  well,  although  assuming  it  consisted  of  copper  ions  with  charge  state  Z,  one 
can  compute  N  ~  1017  from  the  measured  Rb  ^  1cm  (B=0.8T).  Similarly,  one  infers  that 
Pi  =  0.8Z  cm.  In  these  experiments  flute  modes  were  observed  photographically  on  the 
expanding  surface  of  the  plasma  and  the  position  of  the  tips  of  the  flutes  and  the  edge  of  the 
main  plasma  column  could  be  determined  as  a  function  of  time.  The  results  are  shown  in 
Figure  2,  for  various  values  of  the  applied  magnetic  field.  Generally,  one  sees  that  the  tips 
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Figure  2.  [top]  Results  from  Dickinson  ct  al.  [1962]  showing  position  of  flute  tips  (circles) 
and  plasma  core  surface  (x’s)  versus  time  for  various  values  of  the  applied  magnetic 
field,  [bottom]  Onset  time  of  instability  (where  flute  and  core  curves  cross)  versus  B. 
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the  tips  freestream  outward,  while  the  central  core  expands  to  some  maximum  radius  RB. 
If  one  plots  the  time  ( tc )  when  the  flutes  and  core  become  distinct  by  extrapolating  the 


separate  curves  back  until  they  cross,  one  finds  (lower  panel  of  Figure  2)  that  tc  ~  B~0  % 
and  Rc  ~  RB.  It  is  also  interesting  to  note  that  other  studies  discussed  in  Dickinson 
et  al.  [1962]  show  that  the  wavelength  of  the  instability,  A  ~  B~a,a  <  1  (consistent 
with  the  driven  lower  hybrid  drift  instability  [Okada  et  al.,  1981;  Winske,  1988])  and 
Rb  ~  B~s  ,6  <  1  (consistent  with  RB  ~  B~ 2/3  in  3-D). 

D.  Japanese  Laser  Experiments 

The  next  set  of  experim  .nts  involve  a  laser  generated  plasma  expanding  against  a 
strong  magnetic  field,  B  ~  0.5-1.3T  [Okada  et  al.,  1981].  A  glass  laser  (~  100  J)  strikes 
an  aluminum  target,  with  the  plasma  expanding  out  to  RB  ~  l-2cm,  with  pt  ~  RB .  The 
flutes  appear  near  Rc  ~  RB,  but  no  studies  of  the  onset  of  the  instability  with  B  were 
reported. 

E.  NRL  Laser  Experiments 

The  NRL  laser  experiments  [Ripin  et  al.,  1987]  are  similar  to  the  Japanese  experiments 
in  overall  form,  although  the  magnetic  field  was  a  little  weaker  (B=0.1-1.0T)  and  the 
plasma  confinement  radius  was  much  larger,  RB  ~  3  — 15cm,  with  comparable  p,  (assuming 
A/+1°).  In  this  case,  however,  the  instability  appears  well  before  the  plasma  reaches  its 
maximum  confinement  radius.  The  panels  at  the  top  of  Figure  3  show  the  position  of  the 
flutes  (open  circles)  and  plasma  core  (solid  circles)  as  a  function  of  time  for  two  different 
magnetic  field  strengths.  The  points  at  different  times,  however,  are  from  different  laser 
shots,  unlike  the  Dickinson  et  al.  [1962]  experiments  shown  earlier.  However,  similar  to 
Dickinson  et  al.,  the  flutes  freestream  outward,  while  the  plasma  core  comes  to  rest  at 
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Figure  3.  [top]  Results  from  Ripin  et  al.  [1987]  showing  position  of  flute  tips  (open 
circles)  and  plasma  core  surface  (closed  circles)  at  different  times  (from  different  shots), 
[bottom]  Amplitude  of  flutes  for  both  cases  (squares=0.lT;  circles=1.0T)  at  different 
times  showing  lack  of  B  dependence  on  onset  time. 
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R  =  Rg.  For  the  weaker  field  case  (B=0.1T),  Rg  is  about  15  cm,  off  scale  in  thp  picture. 
The  bottom  panel  of  Figure  3  plots  the  data  for  the  flutes  for  both  magnetic  field  values 
(B=0.lT  corresponds  to  the  squares).  Interestingly,  the  data  suggests  that  the  onset  of 
the  instability  is  independent  of  the  magnetic  field,  contrary  to  that  displayed  in  Figure  2. 
F.  Summary 

In  all  four  experiments  a  strong  flute  instability  is  observed.  From  the  data  one  can 
raise  two  important  issues.  The  first  is  where  does  the  instability  turn  on.  When  Rg  is 
large  (AMPTE,  NRL)  the  instability  appears  well  before  Rg.  The  other  two  situations 
(Dickinson,  Okada)  see  the  instability  at  Rc  cz  Rg;  in  these  cases,  however,  Rg  is  only  1-2 
cm  and  perhaps  at  these  small  distances  the  diagnostics  cannot  easily  distinguish  between 
Rc  and  Rg.  The  second  issue  is  how  does  the  turn-on  of  the  instability  depend  on  the 
applied  magnetic  field.  One  knows  that  there  must  be  some  dependence,  as  no  instability 
is  observed  without  the  magnetic  field  being  present  [Dickinson  et  a/.,  1962;  Ripin  et  a/., 
1987],  The  results  of  Dickinson  et  al.  [1962]  further  indicate  (Figure  2)  that  the  turn-on 
depends  strongly  on  B,  tc  ~  B~ 0  8 .  However,  it  should  be  kept  in  mind  that  since  Rc  ~  Rg 
in  that  experiment,  the  turn-on  time  may  just  reflect  the  time  to  maximum  expansion, 
tB  ~  Rb  ~  J3-2/3.  On  the  other  hand,  the  experiments  of  Ripin  et  al.  [1987]  show  little, 
if  any,  dependence  on  B.  A  closer  examination  of  the  bottom  panel  of  Figure  3  shows  that 
during  the  growth  stage  ( t  <  60ns)  the  circles  (stronger  field  case)  tend  to  lie  to  the  left 
of  the  squares,  which  suggests  some  dependence  on  the  field  strength.  Again,  because  the 
points  come  fro-  r  different  shots,  it  is  difficult  to  draw  any  definitive  conclusions  from  the 
figure.  We  thus  turn  to  the  simulations  to  look  for  further  clues. 
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3.  Review  of  Simulations 

A  number  of  simulations  using  different  physics  models  have  been  carried  out  which 
show  the  development  of  flute  modes  on  expanding  plasmas.  We  review  the  relevant  results 
of  each,  then  summarize  at  the  end.  Further  discussion  of  each  of  these  calculations  is  found 
in  Akimoto  et  ah  [1988]. 

A.  Hybrid  Simulations  with  me  =  0 

Brecht  and  Thoms  [1988]  used  a  particle  ion  fluid  electron  hybrid  code  to  study  the 
instability  in  2-D  and  3-D.  In  these  calculations  Px/Rb  is  much  less  than  unity.  The 
instability  appears  when  Rc  cs  RB  in  2-D  and  somewhat  faster  in  3-D.  The  shortest 
allowable  modes  in  the  simulation  dominate.  No  scaling  with  B  has  been  reported. 

B.  Hybrid  Simulations  with  me  ^  0 

Sgro  et  ah  [1988]  employed  a  hybrid  model  containing  electron  inertial  effects.  The 
calculations  were  done  in  a  slab  geometry,  whith  the  plasma  thus  expanding  in  one  di¬ 
rection.  In  these  calculations  the  initial  plasma  state  is  characterized  by  Px/Rb  ~  1  and 
Ti  »  Te.  The  instability  starts  to  develop  before  RB,  but  the  growth  of  the  flutes  is  most 
dramatic  when  the  plasma  recompresses.  Scaling  with  ion  mass,  but  not  B,  was  done.  At 
higher  rriijme  ratio  the  instability  develops  later  in  time,  although  it  is  difficult  to  infer 
an  onset  time  scaling  with  ion  mass. 

C.  MHD  Code  with  Hall  Term 

Huba  et  al.  [1987]  show  that  the  instability  also  occurs  in  a  single  fluid  code,  when  the 
Hall  term  is  included.  (Without  it,  the  usual  Rayleigh-Taylor  instability  grows  instead.) 
The  simulations  were  carried  out  in  a  slab  geometry  by  imposing  an  artificial  gravity;  RB 
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is  not  clearly  defined  in  this  case.  The  instability  grows  rapidly,  at  the  shortest  allowable 
wavelength.  No  scaling  studies  have  yet  been  done. 

D.  Electrostatic  Particle  Code 

Galvez  ct  al.  [1988]  consider  a  hoi  plasma  cylinder  ( R .  >  p,),  expanding  thermally 
against  a  magnetic  field  in  the  electrostatic  limit.  Thus,  no  magnetic  cavity  forms  and 
Rb  is  not  a  relevant  parameter.  It  was  noted  in  these  simulations  that  it  takes  some  time 
before  the  instability  grows,  as  a  charge  layer  is  formed  first.  From  the  results  of  the 
simulations  with  various  values  of  B,  it  is  evident  that  there  is  some  (weak)  dependence 
on  the  magnetic  field,  with  the  instability  developing  sooner  when  B  is  stronger,  although 
a  scaling  of  tc  with  B  has  not  been  extracted. 

E.  Electromagnetic  Particle  Simulations 

Winske  (l988j  used  a  2-D  electromagnetic  particle  code  to  show  the  formation  of 
the  diamagnetic  cavity  and  the  growth  of  short  wavelength  flute  modes.  Figure  4  [from 
Winske  (1989)]  displays  an  example  of  one  such  calculation,  showing  ion  density  contours 
at  various  times.  It  is  clear  that  the  instability  develops  at  Rc  <  Rb  and  the  flutes  stay 
relatively  small  until  the  recollapse  of  the  central  core,  when  then  become  lengthened  as 
the  flutes  continue  to  stream  outward.  The  mode  number  also  changes  at  this  time,  as 
discussed  in  more  detail  by  Winske  [1989].  Runs  with  various  mass  ratios  and  magnetic 
field  strengths  were  also  carried  out.  Figure  5  (again  from  [Winske,  1989])  shows  that  the 
observed  wavenumbers  during  the  growth  stage  (squares)  are  consistent  with  theory  (solid 
lines)  for  the  driven  lower  hybrid  drift  instability.  Of  more  direct  interest  here  is  Figure  6, 
which  shows  the  time  when  the  instability  appears  as  a  function  of  applied  magnetic  field 
and  ion  mass,  tc  ~  B_1/2(mt/me)1/4.  Calculations  in  3-D  (Barnes  ct  al.,  in  preparation] 
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WlJft  m  0 


ULHt  —  4 


Uifit  rn  3  VLffi  =  8 


Figure  4.  Results  of  particle  simulations  [Winske,  1989]  displaying  ion  contours  at  various 
times  to  show  the  expansion  of  the  plasma  and  development  of  the  instability  on  the 
surfs  e. 
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simulations  showing  the  onset  time  of  the  instability  versus 


show  similar  behavior,  and  as  with  those  of  Brecht  and  Thomas  [1988],  indicate  that 
the  instability  develops  somewhat  faster  for  comparable  parameters,  because  the  plasma 
deceleration  is  larger  in  3-D. 

F.  Summary 

The  simulations  show  that  the  instability  does  develop,  in  a  variety  of  physics  models 
consistent  with  linear  theory.  In  most  of  the  simulations  the  instability  turns  on  rapidly, 
i.e.,  Rc  <  Rb ■  They  also  demonstrate  that  the  flutes  do  not  become  large  until  the  core 
plasma  slows  enough  to  become  well  separated  from  the  streaming  flute  tips.  And  the 
simulations  suggest  the  onset  of  the  instability  decreases  with  magnetic  field  strength  and 
increases  with  the  ion  to  electron  mass  ratio.  The  one  exception  is  the  hybrid  simulation 
of  Brecht  and  Thomas  [1988].  In  their  case,  however,  because  p,  -C  i?B,  the  instability  is 
very  weak,  compounded  by  the  fact  that  because  of  the  grid  used  in  the  simulations  the 
wavelengths  are  long,  A  ~  c/ui  (ion  plasma  frequency  =  u>,  =  (47rn,e2/mi)1/2).  Thus,  the 
instability  grows  very  slowly  and  only  appears  when  Rc  c;  Rb,  when  the  deceleration  is 
largest. 

However  given  these  results,  it  is  not  clear  whether  the  simulations  can  say  anything 
about  the  actual  turn-on  radius  Rc  or  turn-on  time  tc  in  an  experiment.  This  is  due  in  part 
because  the  simulations  start  with  an  ideal  plasma  state  that  has  a  well  defined  shape  and 
velocity  distribution  that  leads  to  rapid  growth  of  the  instability.  The  initial  conditions  in 
the  actual  experiment  are  much  more  complex,  because  the  plasma  is  generated  by  a  laser 
(or  gun)  beam  striking  a  target,  or  by  a  release  of  m  •  ral  atoms  that  ionize  over  a  finite 
time  interval.  In  addition,  because  of  computational  constraints,  the  initial  radius  of  the 
plasma  R„  is  not  very  small  in  relation  to  Rc  and  Rb ,  rather  R*  ~  Rc  <  Rb .  A  third 
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complication  is  the  presence  of  a  background  plasma  in  the  actual  experiments,  absent  in 
most  of  the  simulations.  In  the  laser  experiments,  for  example,  the  ionization  level  of  the 
background  may  be  a  strong  function  of  position  [Brecht,  private  communication,  1988], 
which  may  affect  the  turn-on  position. 


4.  Theoretical  Scalings 
A.  Preliminaries 

We  consider  the  driven  lower  hybrid  drift  instability  in  the  limit  of  cold  ions,  in  which 
case  the  relative  electron  ion  streaming  is  due  solely  to  the  deceleration  of  the  ions  by  the 
magnetic  field  [Winske,  1988].  The  deceleration  g  is  defined  as 
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dVXT 

dt 


RbKRb} 


(5) 


where  a  =  1,5  =  1(2-D)  or  a  =  2 ,6  =  3/2  (3-D);  g  thus  increases  with  radius  and  is  a 
maximum  at  R  —  Rb.  The  above  expression  assumes  an  expansion  perpendicular  to  B; 
in  3-D  one  can  also  include  the  free  expansion  along  the  magnetic  field,  which  changes 
6  — »  2.3  [Gisler  and  Lemons,  1988], 

The  corresponding  cross-field  drift  is 


V,  =  5/n,  (6) 

and  the  actual  quantity  of  interest  is  V8/va  ,  where  va  is  the  Alfven  speed,  =  jB2/47rnm, 
(va  is  a  function  of  R,  because  the  density  n  is). 

We  assume  two  different  models  for  n(R).  First,  the  plasma  expands  as  a  uniform 
cylinder: 

N  =  =  n{R)R2Ln  (7) 


IB 


where  n, ,  Rt  are  the  initial  density  and  radius  respectively.  Second,  the  plasma  compresses 
into  a  thin  shell  A,  such  that  n(R)  =constant,  i.e., 


N  =  n(R)2nRAL  =>  A  ~  1/7? 


The  first  model  implies 


Vg/vA  —  Pi/Rb  —  constant  (2  —  D) 


(8) 


(9) 


and 


^  =  (3-D> 


while  the  second  yields  (in  2-D) 


Vg/vA  =  (R/'RbKpx/Rb) 


(10) 


(11) 


i.e.,  Vg/vA  increases  linearly  with  R.  In  3-D,  this  same  model  gives 

v,/»A  =  IwRbVMRb)  (12) 

Results  of  2-D  particle  simulations  [Winske,  1988]  in  which  Vg/vA  is  measured  at  various 
R  shows  behavior  consistent  with  (11)  (Figure  7). 

B.  Onset  Criteria 

We  consider  four  different  criteria  for  the  turn-on  of  the  instability,  deriving  each  for 
one  particular  model  [2-D  expansion,  n(R)=constant].  The  following  subsection  summa¬ 
rizes  the  results  for  all  the  models. 

The  first  criterion  is  that  the  instability  turns  on  when 


Vg/vA  >  2 


(13) 
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When  the  ions  are  hot  (ft  ~  1  =»  vA  ~  vt),  (13)  corresponds  to  Vg  >  v,,  the  nonresonant 
regime  of  the  lower  hybrid  drift  instability.  The  condition  (13)  is  also  analogous  to  the  one 
derived  by  Huba  [1987]  for  the  unmagnetized  ion  Rayleigh  Taylor  instability.  In  that  case, 
the  assumption  of  unmagnetized  ions  implies  7  >  fl;.  Hassam  and  Huba  [1987]  show  that 

7  =  k(g/ en)1/2  (14) 


where  k  is  the  wavenumber  and  en  =  n~1dnfdx  is  the  density  inverse  scale  length.  Because 
g  ~  R g3  ~  B2,  the  righthand  side  of  (14)  is  proportional  to  B  and  the  expression  7  >  f 1, 
is  independent  of  B. 

Our  criterion  is  based  on  the  results  of  Figure  8,  showing  the  growth  rate  maximized 
over  k  and  the  corresponding  wavenumber  versus  Vg/vA  ( Te  =  0 ,ft  =  10_4,enc/u;,  =  10). 
For  Vg/vA  <  2,  both  7  and  k  change  rapidly,  while  for  Vgiv\  >  2,  they  vary  much  more 
slowly.  The  condition  (13)  can  thus  be  written  as 


Yl 


R  pi 

R-b  Rb 


>  2 


(15) 


or 


Rc  >  2 Rb/{p%/Rb) 


(16) 


Normalizing  to  the  reference  case  where  B  =  B0,  Rb  =  R0,  =  m0  and  using  the  fact 

Rb/Ro  ~  (B0/ B) (m,/m0) 1/2,  pijp0  =  (B0/ B)(mi/m0)  [in  2-D,  from  Eq.  (1-3)] 


Rc  ^  Bp  2 

Ro  B  ( p0 / Ro ) 


(17) 
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Figure  8.  Results  of  linear  theory  for  the  driven  lower  hybrid  drift  instability  sh 


growth  rate  maxiznimized  over  wavenumber  and  the  corresponding  k  vers 


(fit  =  lo-VncM  =  10). 
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The  second  criterion  is  based  on  the  linear  growth  rate.  The  waves  grow  from  noise 
starting  at  R  =  R\  and  are  seen  when  they  reach  a  given  size 

h  =  height  of  flute  =  She11  (18) 


where  6h  is  the  initial  perturbation  (noise).  But  t  ~  R/Vd  and  7  ~  ulh  [lower  hybrid 
frequency  =  ulh  —  (m^ / me ) 1  /2 [  so 


Rc  =  —In  (h/Sh) 
ulh 


(19) 


or 


Rc  B0  /  rn,  \  1/2 

R0  B  m0 

The  third  turn-on  condition  is  based  on  the  fact  that  because  the  initial  density 
cloud  is  so  large  the  electron  thermal  beta,  /?e  =  87 meTe/'B2,  is  much  greater  than 
Because  the  driven  lower  hybrid  drift  instability  is  stabilized  by  0t  effects  [Drake 
1983],  the  instability  turns  on  when  0e  drops  to  some  low  value,  ~  1: 


(20) 

of  the 
unity. 
et  ai, 


0, 


8nneTe 

B2 


(21) 


If  ne  =constant,  because  of  compression  of  the  cloud  into  a  thin  shell,  (3e  ^constant  and 
this  condition  does  not  apply.  On  the  other  hand,  if  we  use  the  other  model,  n(i?)  ~  R~2 
(21)  becomes 

where  0O  =  0e{Ro). 

The  last  criterion  for  instability  onset  is  again  based  on  Figure  8,  which  shows  that 
for  large  Vg/v a,  k  does  not  change  very  rapidly.  As  the  plasma  radius  R  increases,  Vg/vA 
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increases  so  that  k  decreases.  It  is  assumed  that  the  instability  will  grow  when  the  mode 


number  (m=kR)  changes  slowly  in  one  growth  time.  Figure  8  (and  the  linear  analysis  in 
[Winske,  1988])  indicate  that  for  Vg jvA  >  1, 

~  (Vs/va)-113  (23) 


so  that 


using  (11).  Thus 


m  =  kR 


Rn(R)1/2 

(va/vAy/* 


A  m/m  ~  A  RjR 


(24) 


(25) 


and 


A R  ~  V{R)At  s  VD{\- R2  jRlY^At 


(26) 


Taking  At  ~  -y  1  ~  we  thus  find 


Am  ^  Pi  (1  -  R2/R2B)l/2  me  1/2 
m  ~  Rb  R/Rb  1  q 


m. 


(27) 


where  q  is  some  small  number,  e.g.,  0.1.  Then  one  can  easily  show  that  R/Rb  is  small  so 
that  (1  -  R2 /R2b)x^2  can  be  neglected,  yielding 


Re 


lc  I  /  me  .  i /2 ,  m,  ^ yg  Po  Bp 

R0  q  rn0  m0  5 


(28) 


C.  Results  for  All  Models 

The  results  for  both  density  models  [Model  1:  n(R)  ~  R~a  (a  =  2  (2-D)  or  3(3-D)) 
or  Model  2:  n(R)  ~constant]  in  both  2-D  and  3-D  are  summarized  in  the  following  table. 
We  express  the  scaling  with  magnetic  field  and  ion  mass  in  the  following  form: 


R 

Rc 


(frer 


(29) 
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2-D  3-D 


Criterion  # 

Model  1 

Model  2 

Model  1 

Model  2 

1 

— 

0=  -1 

0  =  — 1/3 

0  =  — 1/2 

p  =  0 

P  =  —1/3 

p  =  0 

2 

0  =  ~1 

0  =  -l 

/?  =  -! 

/?  =  -! 

M  =  1/2 

. 

P  =  l/2 

P=  1/2 

P  =  1/2 

3 

0=~1 

— 

0  =  -2/3 

— 

p  =  0 

t 

p  =  0 

4 

— 

0  =  ~1 

0  =  ~l 

1 

0  =  -l 

P=  1/2 

P  =  1/2 

P  =  1/2 

D.  Interpretation 

We  interpret  the  results  displayed  in  the  above  table  as  follows.  Conditions  2  and 
4  essentially  depend  on  the  linear  growth  rate,  7  ~  win,  and  not  on  the  density  model 
assumed.  So  for  all  four  models  when  these  two  conditions  apply,  one  has  R/Rc  ~  ~ 

(B / Bo)"1  (mi/rrio)1/2 .  The  third  condition  depends  on  0e,  hence  has  no  dependence  on 
ion  mass;  when  the  density  does  not  vary  with  radius,  neither  does  0e,  and  this  condition 
is  not  applicable.  Only  the  first  turn-on  condition  shows  significant  model  variation  of  the 
B  and  m,  dependence,  and  displays  the  weakest  B  dependence  in  3-D.  None  of  the  2-D 
conditions  show  the  same  scaling  as  seen  in  the  2-D  simulations  [Figure  6].  The  fact  that 
the  simulation  scaling  is  consistent  with  with  a  ~  1/2,  could  suggest  that  conditions 
2  or  4  may  be  applicable,  but  that  the  asymptotic  value  of  the  growth  rate  (7  =  u>lh) 
is  not  quite  correct.  None  of  the  criteria  derived  here  have  no  dependence  on  the  applied 
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magnetic  field,  as  observed  by  Ripin  et  al.  [1987],  although  the  dependence  in  one  case  is 
fairly  weak  (Condition  1,  Model  1  in  3-D).  We  will  discuss  these  results  further  in  relation 
to  the  experiments  next. 

5.  Discussion 

Previously,  we  have  argued  that  the  fluting  which  appears  on  the  surface  of  plasma 
clouds  expanding  across  an  ambient  magnetic  field  is  due  to  the  driven  lower  hybrid  drift 
instability  [Winske,  1988].  Because  this  is  a  high  frequency  short  wavelength  mode,  to 
explain  some  of  the  observable  manifestations  of  the  instability,  nonlinear  effects  need  to 
be  invoked  [Winske,  1989].  In  this  report  we  have  considered  conditions  for  the  turn-on  of 
the  instability  deduced  from  linear  theory  and  apply  them  to  various  numerical  simulations 
and  experiments  in  which  the  instability  is  excited.  Four  different  onset  conditions  have 
been  derived,  and  each  has  been  worked  out  for  both  2-D  and  3-D  expansions  assuming  two 
different  models  for  the  plasma  density  profile.  This  range  of  models  has  been  considered 
in  order  to  see  if  there  is  much  sensitivity  of  the  magnetic  field  and  ion  mass  dependence 
on  the  model  used.  In  addition,  the  various  simulations  and  experiments  tend  to  lie  in 
different  parameter  regimes,  and  thus  we  do  not  expect  all  the  data  to  be  explained  by  a 
single  scaling  law. 

Generally,  the  simulations  are  easier  to  understand.  This  is  because  the  initial  con¬ 
ditions  are  well  known  and  one  can  (in  principle)  determine  local  quantities,  such  as  drift 
velocities,  density  gradients,  temperatures,  etc.,  as  needed  by  the  theory.  In  addition,  the 
simulations  are  free  of  complications  such  as  changing  charge  states,  impurities,  radiation, 
background  ionization  effects,  etc.  The  simulations  in  which  scaling  of  the  instability  onset 
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has  been  attempted  show  faster  turn-on  with  stronger  magnetic  field  and  smaller  ion  to 
electron  mass  ratio.  These  same  trends  are  predicted  for  the  linear  growth  rate  of  the 
instability.  Thus,  the  simplest  consistent  explanation  is  that  the  onset  time  tc  scales  as  the 
inverse  growth  rate,  7-1.  That  the  simulations  do  not  follow  a  simple  scaling  like  (29)  with 
j3  =  —  1  and  p  =  1/2  just  reflects  the  fact  that  the  growth  rate  used  in  those  expressions 
(7  ~  ojlh)  is  the  idealized  limit.  Of  course,  it  should  also  be  kept  in  mind  that  the  initial 
conditions  in  the  simulation  are  chosen  to  promote  rapid  growth  of  the  instability. 

The  experiments,  however,  are  more  ambiguous.  From  the  published  data  of  the 
plasma  gun  experiment  of  Dickinson  et  al.  [1962]  one  can  derive  an  onset  time  that  scales 
inversely  with  B,  tc  ~  B~0  8.  One  possible  explanation  is  that  in  this  experiment  the 
instability  is  strongly  driven,  given  the  large  flutes  that  are  observed,  i.e.,  Pi/ Rb  >  1-  The 
strong  dependence  on  B  then  reflects  the  large  linear  growth  rate,  7  ~  ulh-  But  assuming 
a  copper  plasma  and  B=0.8T,  one  readily  computes  ~  2ns  (assuming  Z=l).  Because 
the  experimental  onset  time  at  B=0.8T  is  about  1.2psec,  this  implies  about  103  e-folding 
times.  Obviously,  there  is  plenty  of  time  for  short  wavelength  modes  to  grow,  saturate, 
and  nonlinearly  evolve  into  the  structures  actually  seen.  Alternatively,  some  mechanism 
such  as  finite  /?e  effects  prevents  the  instability  from  growing  until  much  later,  t  ~  1/xsec. 
Another  possibility,  already  mentioned,  is  that  the  tc  scaling  reflects  the  fact  that  the 
instability  onset  occurs  near  Rb,  hence  tc  ~  ts  ~  Rb  ~  B~2^3. 

In  the  AMPTE  experiment,  on  the  other  hand,  one  does  not  have  the  problem  of 
trying  to  decipher  a  25  year  old  experiment,  but  as  there  were  only  two  magnetotail 
barium  releases,  one  does  not  have  a  scaling  for  the  onset  of  the  instability.  However,  one 
is  again  faced  with  the  fact  that  the  onset  time  ~  180sec  is  much  greater  than  ~  2 
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sec,  suggesting  the  instability  is  inhibited  from  growing  at  early  times.  This  could  be  due 
to  the  fact  that  the  plasma  beta  is  initially  very  large  and  the  background  plasma  is  hot,  so 
that  finite  (3e  effects  can  be  used  to  explain  the  slow  onset.  Or  again,  because  of  the  large 
discrepancy  between  the  observed  waves  and  those  predicted  for  the  driven  lower  hybrid 
drift  instability,  it  is  conceivable  that  the  short  wavelength  modes  have-  grown  up  and 
saturated  and  one  is  seeing  the  slower  growing,  nonlinear  manifestation  of  the  instability 
(in  which  case  the  onset  criteria  derived  here  are  probably  not  relevant). 

That  leaves  the  laser  experiments,  especially  the  experiments  of  Ripin  ef  al.  [1987 j 
in  which  the  onset  time  was  found  to  be  independent  of  B.  For  Al  +  1°  in  a  0.1T  magnetic 
field,  ore  finds  =  4ns  compared  to  tc  =  50ns.  There  are  many  e-foldings  possible 
(and  many  more  for  the  1.0T  field  case),  again  suggesting  some  initial  inhibition  of  the 
instability  (which  could  be  finite  or  complications  due  to  the  changing  charge  state 
of  the  aluminum  or  the  background  plasma  conditions).  Because  the  experiment  is  still 
operational,  it  suggests  some  further  avenues  of  investigation  to  help  clarify  the  onset 
question.  With  the  improved  capability  to  obtain  multiple  exposures  per  shot,  it  would  be 
most  useful  to  again  check  the  B  dependence  of  the  onset  time,  to  see  if  some,  even  weak, 
dependence  on  the  magnetic  field  exists.  Some  of  the  theoretical  onset  conditions  predict 
weaker  magnetic  field  dependence  in  a  realistic  three  dimensional  expansion  compared  to 
the  two  dimensional  simulations.  A  second  set  of  useful  experiments  would  be  to  vary  the 
target  material  to  test  the  mass  scaling.  No  dependence  on  the  ion  mass  could  indicate  that 
finite  electron  beta  effects  inhibit  the  instability  growth  at  early  times.  A  weak  dependence 
on  mass,  particularly  if  opposite  to  a  slightly  stronger  dependence  on  B,  could  suggest  a 
growth  rate  type  of  scaling  similar  to  that  seen  in  the  simulations.  A  final  useful  addition 
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to  the  experiment  would  be  a  more  thorough  study  of  the  data  from  the  electrostatic  probe 
that  detects  the  growth  of  the  very  short  wavelength  high  frequency  noise.  Scaling  of  that 
noise  with  B  and  could  show  that  lower  hybrid  modes  are  indeed  present  and  suggest 
that  the  actually  observed  waves  are  a  nonlinear  evolution  of  these  shorter  wavelength 
modes. 

As  a  final  note,  it  should  be  mentioned  that  not  all  of  the  burden  falls  on  the  ex¬ 
perimenters.  The  simulators  should  also  direct  their  efforts  toward  making  the  numerical 
calculations  more  realistic,  i.e.,  including  3-D  effects,  improving  the  initial  conditions  of 
the  plasma,  and  modeling  some  of  the  messier  details  (charge  states,  incompletely  ionized 
background)  of  the  laser  experiment  to  assess  how  important  a  role  they  play. 
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